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AN  EXPERIMENTAL  DETERMINATION  OF  THE  TRANSITION  TO 
TURBULENCE  USING  THE  ELECTRIC  DOUBLE  LAYER 
By 

August,  1965 

Chairman:  K.  T.  Millsaps 

Major  Department:  Aerospace  Engineering 

The  experimental  determination  of  the  critical 
Reynolds  number  for  transition  from  laminar  to  turbulent 
flow  in  a straight  circular  pipe  has  been  accomplished  by 
use  of  a streaming  potential  technique.  The  streaming 
potential  was  generated  by  an  electrolyte  (0.001  M KC1) 
flowing  through  a capillary  tube  with  a diameter  of  0.0242 
Inches.  It  was  found  that,  at  the  onset  of  turbulence, 
the  streaming  potential  began  to  fluctuate  in  an  easily 
identifiable  manner.  This  fluctuation,  combined  with  mass 
flow  measurements,  was  used  to  determine  the  critical 
Reynolds  number  to  be  1907  ± 3 for  a fixed  geometry  and 


INTRODUCTION 


The  purpose  of  the  present  Investigation  is 
twofold.  First,  a new  technique  using  the  "electric 
double  layer"  is  developed  which  is  capable  of  measur- 
ing hydrodynamic  effects  within  Angstroms  of  a wall. 
Secondly,  this  method  is  applied  to  the  old,  but  as 
yet  theoretically  unresolved,  problem  of  the  transition 
from  laminar  to  turbulent  flow  in  a straight  pipe  of 
circular  cross-section. 


2.  Stability  of  Poiseullle  Flow 

The  theoretical  problem  of  determining  the 
stability  of  Poiseullle  flow  has  been  the  subject  of 
investigations  for  many  years;  considerable  efforts 
have  been  made  for  both  the  plane  and  circular  cases. 
Interesting  results  have  been  obtained  in  both  cases 
although  the  former  case  has  been  by  far  the  most 
fruitful. 

Plane  Polseuille  flow  has  been  studied  on 
the  basis  of  two  dimensional  infinitesimal  disturbances 


where  the  perturbation  in  the  stream  Junction  is  of 

♦ fcc.Y.t)  =4>(r)e"“ft"ct) 

where  ')/'  is  the  perturbation  in  the  stream  function, 

tj>  is  the  amplitude, 
cc  is  the  wave  number, 

C is  the  speed  of  propagation. 


The  rectangular  Cartesian  coordinates 


parallel  and  perpendicular  to  the  flow.  This  form  is 
not  quite  as  restrictive  as  it  might  seem  since  Squire 
has  shown  that  the  three  dimensional  case  is  equivalent 
to  the  two  dimensional  one  at  a lower  Reynolds  number. 
Heisenberg*2'  was  the  first  to  deduce  the  instability 
of  plane  flott  for  sufficiently  large  Reynolds  number , 
but  he  did  not  calculate  a critical  value  beyond  which 
instability  would  occur.  The  method  of  Heisenberg 
was  extended  by  Tollmien*5'  and  Schlichting  and  was 
used  to  calculate  the  stability  curves  for  the  boundary 
layer  on  a flat  plate.  Later,  Lin15'  clarified  the 
previous  work  and  succeeded  in  calculating  the  neutral 


plane  Poiseuille 


Contrary 


conclusion  of  Lin,  Pekeris*6*  applied  a different 
technique  and  concluded  that  Poiseuille  flow  was 
stable  for  all  Reynolds  numbers.  To  resolve  this 
disagreement,  Thomas*7*  calculated  the  critical 
Reynolds  number  by  direct  numerical  methods  and  found 
it  to  be  5780.  However,  one  should  bear  in  mind  that 


the  above  work  is  for  infinitesimal  two  dimensional 
disturbances  and  that  nothing  has  been  successfully 
worked  out  for  an  arbitrary  perturbation. 

The  results  for  pipe  flow  are  not  as  rewarding 
as  the  plane  case.  Sexl*8)  was  the  first  to  solve  the 
viscous  stability  equations  for  axisymmetric  disturbances. 
However,  his  results  are  questioned  because  of  the  math- 
ematical simplifications  used.  Pretsch*8*  found  that  if 
the  analysis  is  restricted  to  a thin  layer  near  the  wall 
the  problem  becomes  the  same  as  for  a disturbance  near 
the  wall  in  Couette  flow.  Pekeris*10*  found  a different 
solution  for  a region  near  the  center  of  the  pipe.  Corcos 
and  Sellars*  * have  obtained  a solution  which  accounts 
for  the  work  of  both  Pretsch  and  Pekeris;  however,  their 
solution  is  in  terms  of  a finite  set  of  eigenfunctions. 
Schensted*12’  has  resolved  this  question  and  has  shown 
that  for  axisymmetric  disturbances  the  set  of  eigenfunctions 


are  not  only  Infinite  but  also  complete.  All  the  existing 
investigations  of  the  parabolic  profile  tend  to  show  that 
the  flow  is  stable  with  respect  to  infinitesimal  two 
dimensional  disturbances.  Of  course,  it  is  well  known 
that  turbulent  flow  does  exist  in  an  experimental  sense  in 
a circular  pipe.  Two  possible  explanations  of  this  paradox 
have  been  proposed;  one  by  Meksyn  and  Stuart* 13^  and  the 
second  by  Tatsumi**3 4^ . Meksyn  and  Stuart  showed  that  in 


a channel  the  minimum  critical  Reynolds  number  decreased 
as  the  amplitude  of  the  oscillations  increased,  thus 
suggesting  that  instability  in  a circular  pipe  is  the 
result  of  finite  oscillations.  On  the  other  hand,  Tatsumi 
considered  the  stability  in  the  inlet  section  of  the  pipe 


where  the  flow  is  more  of  the  boundary  layer  type;  he 
calculated  a minimum  critical  Reynolds  number  of  the 
order  of  10  , whereas  the  maximum  Reynolds  number  which 
has  been  obtained  experimentally  by  Taylor  (unpublished) 
and  Ekman  for  laminar  flow  is  4 x 10  . 


3-  Experimental  Investigations  of  Transition  in  a Pipe 
The  experimental  study  of  the  stability  of  flow 

through  pipes  of  circular  cross-section  has  occupied  a 
prominent  position  in  the  history  of  fluid  mechanics 


since  the  time  of  Osborne  Reynolds.  In  1883,  Reynolds' 
reported  the  results  of  experiments  which  showed  beyond 
doubt  that  there  exist  two  possible  modes  of  flow  in  the 
same  pipe — laminar  and  turbulent.  In  addition  to  the 
experiments,  Reynolds  reasoned  dimensionally  that  the 
flow  in  pipes  of  different  diameters  should  be  similar 
if  the  dimensionless  number  (Reynolds  number)  Z^4 ^ 
is  the  same  in  each  pipe  where  yO  is  the  density  of 


a characteristic 
an  velocity  and  u. 


3 independent  s 


linear  dimension,  it  is  the  me 
is  the  absolute  viscosity. 

Reynolds  performed  t 
experiments.  In  the  first  group,  which  is  the  most 
widely  known,  he  Introduced  dye  into  the  mouth  of  a 
tube  while  the  velocity  was  increased.  He  found  that 
at  low  velocities  the  dye  would  go  through  the  tube  as 

the  velocity  increased,  a point  would 
e filament  would  very  suddenly  spread 


a filament,  but, 
be  reached  where 
out  and  fill  the  downstream  end  c 
he  called  the  critical  velocity, 
measurements  are  best  summarized 
"Those  experiments  had 
that  there  was  a critical  va 
at  which  eddies  began  if  the 
mately  steady  when  drawn  lnt 


rhe  results  of  these 
r Reynolds: 
lown  definitely 
le  of  the  velocity 
rater  were  approxi- 


had  also  shown  definitely  that  at  such  critical 
velocity  the  water  in  the  tube  was  in  a highly 
unstable  condition,  any  considerable  disturbance 
in  the  water  causing  the  break  down  to  occur  at 
velocities  much  below  the  highest  that  could  be 
attained  when  the  water  was  at  its  steadiest; 
suggesting  that  if  there  were  a critical  velocity 
at  which,  for  any  disturbance  whatever,  the  water 
became  stable,  this  velocity  was  much  less  than 
that  at  which  it  would  become  unstable  for 
infinitely  small  disturbances;  or,  in  other  words 
suggesting  that  there  were  two  critical  values 
for  the  velocity  in  the  tube,  the  one  at  which 
steady  motion  changed  into  eddies,  the  other  at 
which  eddies  changed  into  steady  motion." 

Reynolds  proceeded  to  determine  this  lower  critical 
velocity.  The  experiment  was  to  introduce  initially 
turbulent  water  "from  the  Manchester  main"  into  a long 
straight  lead  pipe  and  to  measure  the  pressure  loss  over 
a 5 foot  length  on  the  downstream  end.  The  critical 
velocity  was  taken  as  the  point  where  the  pressure 
gradient  was  no  longer  proportional  to  the  velocity.  His 
results  showed,  in  modern  terminology,  a lower  critical 
Reynolds  number  of  approximately  2000. 

Since  1883,  many  workers  have  repeated  Reynolds' 
experiments;  the  most  extensive  repetition  of  Reynolds' 


work  was  conducted  by  Stanton  and  Pannell<17\  In  addition, 
several  Interesting  variations  of  the  classical  method 
have  been  devised.  For  example,  Barnes  and  Coker*18^ 
used  a thermal  scheme  where  the  walls  of  the  tube  were 
heated  and  the  onset  of  turbulence  was  detected  by 
noting  a sharp  rise  in  the  temperature  in  the  interior 
of  the  flow.  They  found  the  critical  Reynolds  number 
equal  to  1900.  Also,  in  1962,  Reiss  and  Hanratty<19) <20> 
developed  a technique  from  which  they  could  infer  the 
behavior  of  the  flow  within  the  so-called  "laminar 
sublayer"  by  measuring  the  mass  transfer  to  a small  sink 
on  the  wall.  They  used  a polarized  electrode,  which  was 
current  limited  by  mass  transfer,  as  a sink.  Thus,  any 
fluctuation  in  the  concentration  would  result  in 
fluctuations  at  the  wall.  Using  this  scheme,  they  found 
random  fluctuations  within  the  sublayer  and  from  the 
onset  of  these  fluctuations  deduced  the  critical  Reynolds 
number  to  be  2100. 


4.  The  Electric  Double  Layer 

It  is  well  known  from  the  study  of  colloids 
that  at  a solid-liquid  interface,  the  solid  surface 
becomes  charged.  The  charges  on  the  wall  are  supplied 


liquid 


attached 


where  they, are  attached  is  still  open  to  question. 
Therefore,  following  Helmholtz^21^ , the  experimental 
evidence  that  such  a charged  layer  exists  shall  be 
accepted  in  the  present  work.  Since  there  is  a charge 
residing  on  the  interface,  there  must  also  be  an  equal 
and  opposite  layer  of  charge  in  the  adjacent  region. 
These  charged  layers  are  called  the  electric  double 

The  early  workers^ 21  ^22^  in  the  field 
considered  the  double  layer  to  be  composed  of  two 
distinct  layers  of  charge;  one  fixed  to  the  wall  and 
the  other  free  to  move  with  the  fluid.  Gouy<23)  and 
Chapman*24^  independently  formulated  the  theory  of  the 
diffuse  double  layer  where  the  charge  in  the  fluid  is 
a function  of  the  distance  from  the  wall.  This  is  in 
essence  the  theory  of  ionic  atmospheres  given  some  ten 
years  later  by  Debye  and  Hfickel*25* . Gouy  and  Chapman 
considered  the  double  layer  to  be  ions  attached  to  the 
solid  phase  and  an  equivalent  ionic  charge  of  opposite 
sign  distributed  in  the  solution  near  the  interface. 

uniform  surface  charge.  The  space  charge  in  the  solution 


is  due  to  the  unequal  distribution  of  ions  which  are 
assumed  to  have  no  physical  size.  The  solvent  is 
considered  as  a continuous  medium  entering  the  problem 
only  through  its  dielectric  constant.  Coulomb  inter- 
actions in  the  system  are  described  by  Poisson's 
equation 

v*y  = .iZZ  , 

where  Y'  is  the  potential  (having  a value  of  V»  at 
the  wall) , 

Y is  the  charge  density, 

£ is  the  dielectric  constant, 
and  71  is  the  Laplaclan  operator. 

The  number  density  of  the  ith  ion  species  is  assumed 
to  be  given  by 


r>i  =m.  cxpl-iieY/kT) 

D;  is  the  number  density  of  the 


is  the  valence, 
e is  the  electronic  charge. 


k is  Boltzmann's  constant, 
and  T is  the  temperature  of  the  solution. 

The  charge  density  is  given  by 

7 ' = Z^en; 

Combining  equations  (1),  (2),  and  (3),  one  obtains 

vV  = -f?2;en ;0  exp  fa;  cV/lcT) 

which  is  the  differential  equation  for  the  potential  as 
a function  of  the  space  coordinates.  If  the  assumption 
of  an  infinite  flat  wall  is  made,  equation  (4)  reduces 

expfac  VAt) 

which  is  subject  to  the  boundary  conditions 


distance 


Letting  7C*=  y‘^n'*i‘ 

the  solution  of  (5)  for  a single  species  is 

-y  - | [exp(*c¥'/Zkr)  + l][oipfeeV/kr)-l1 
[exp  (let/ikT)  - ( ] IcxpCzeY./fcT)  *l] 


This  can  be  simplified,  if  ief/kT 
by  expanding  the  exponentials.  The 

Xx  = lnryr 


is  assumed  small, 


f=f.e 


shoving  that  3^"  drops  by  a factor  of  e~'  in  the 
distance  '/*  . This  value  is  called  the  thickness 

of  the  double  layer  or  Debye  length.  As  an  example, 
this  distance  is  approximately  10~  cm.  for  a 0.001  M, 

KC1  solution.  This  theory  vas  found  to  be  an  over 
simplification  since  the  concentration  of  the  electrolyte 
is  calculated  to  be  300M  at  the  vail  for  a 0.1  M solution 
and  a vail  potential  of  200  mv.  This  is  an  impossible 


The  Gouy-Chapman  theory  was  modified  by 
Stern*26^  to  allow  for  the  finite  size  of  the  ions  and 
thus  to  avoid  the  difficulty  mentioned  above.  Stern 
suggested  that,  since  the  influence  of  the  dimensions 
of  the  ions  is  greatest  at  the  wall,  it  is  perhaps 
possible  to  add  a major  correction  to  double  layer 
theory  by  accounting  for  the  finite  size  of  the  ions 
only  in  the  region  of  the  first  few  molecular  layers. 
This  simply  requires  that  the  Gouy-Chapman  theory  be 
applied  not  at  the  wall  but  a short  distance,  & , 

away  from  it.  In  addition,  Stern  considered  the 
possibility  of  adsorption  of  the  ions  in  the  plane 
X - £ . This  layer  has  become  known  as  the  Stern 

5.  Streaming  Potential 

The  streaming  potential  is  caused  by  the  flow 
of  an  electr.olyte  sweeping  out  ions  in  the  movable  part 
of  the  double  layer.  This  flow  causes  an  electrical 
current  by  the  convection  of  charges  parallel  to  the  wall 
which,  in  turn,  results  in  a potential  difference  between 
the  two  ends  of  the  passaged  ' This  potential  difference 
causes  a conduction  current  to  flow  which  at  equilibrium 
is  equal  and  opposite  to  the  convection  current. 


Consider  n tube  of  arbitrary  cross-section  with 
coordinates  as  shown  In  Fig.  1.  Assuming  the  charge 
density  on  the  wall  is  Independent  of  the  flow  velocity, 
the  convection  current , Xt  , one  finds  that 


= f^ir -a  dA 


here  7*  is  the  charge  density, 

V is  the  hydrodynamic  velocity, 

H is  the  unit  normal  to  dA , 

A is  the  cross-sectional  area, 
d R is  equal  to  dNds , 

N is  the  coordinate  normal  to  the  wall , 

f Poisson's  equation  is  used,  equation  (8)  becomes 


d,  =-fjf!^VvndA 


It  is  now  assumed  that  the  thickness  o 
is  small  compared  to  the  dimensions  of 
this  case  V becomes  ,7/p*- 


e double  layer 


= ~$f(  J&VJA 


Fig.  1 Coordinate  system  for  tube 
of  arbitrary  cross-section. 


ipplic 


N-0 


,,  = dHdS 

= -frrfv?  Jh" ds  ~iifjjfirfidHds 


(U) 


...r.  t-rwL. 


dKI  - X, 
JN  L:o  '/T 


Therefore  equation  (11)  becomes 

da] 

The  second  term  in  equation  (12)  will  be  neglected  since 
it  may  be  shown  that  the  ratio  of  the  second  term  to  the 
first  is* 

o(^) 

pressure  gradient  is  assumed  constant  over  the  cross- 
section,  equation  (12)  may  be  written  as 


= t_K 


The  conduction  current,  Xi  , is 


^ is  the  specific  conductance 
E is  the  potential  caused  by  the  fluid  flow. 


The  notation  OM 


sectional  area,  (14)  becomes 


Equilibrium  requires  ^ = -Xi  ; hence 

d£  _ 6* 

dP  4/TjH  -l 

Under  the  above  assumptions,  equation  (16)  shows  that 
the  streaming  potential  is  a linear  function  of  the 
pressure  drop  along  the  tube.  The  first  term  on  the 
RHS  of  equation  (11)  indicates  that  the  convection 
current,  and  hence  the  streaming  potential,  should  be 
sensitive  to  the  velocity  gradient  at  the  wall.  Thus 
a change  in  the  streaming  potential  should  provide  a 
means  of  detecting  a change  in 
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II. 


EXPERIMENT 


1.  General  Description 

The  theory  given  in  section  5 of  Chapter  I 
shows  that  the  streaming  potential  varies  linearly  with 
the  pressure  gradient  along  the  passage.  Also, 
equation  (11)  shows  that  the  streaming  potential  varies 
linearly  with  the  velocity  gradient  at  the  wall.  It  is 
this  relation  that  has  produced  a new  technique  for 
hydrodynamic  measurements . 

Let  us  now  restrict  the  discussion  to  the  flow 
in  a straight  tube  with  a circular  cross-section.  If 
the  flow  is  laminar,  the  velocity  distribution  becomes 
U =4/Ldi  (/-&) 

at  a sufficient  distance  from  the  inlet;  where 

p is  the  pressure, 

^ is  the  distance  along  the  tube, 

/l  is  the  distance  from  the  center  of  the  tube, 


velocity. 


If  this  relation  is  substituted  into  equation  (II),  the 
streaming  potential  is  again  found  to  vary  linearly  witl 
the  pressure  gradient  along  the  tube.  This  has  been 
observed  by  physical  chemists  for  many  years'28^.  If  tl 
flow  changes  from  laminar  to  turbulent,  the  velocity 
gradient  at  the  wall  no  longer  varies  linearly  with  the 
pressure  gradient,  and  one  therefore  expects  a change 
from  a straight  line  to  a nonlinear  curve  when  plotting 
streaming  potential  versus  pressure  drop  along  the  tube 
In  addition,  the  streaming  potential  might  be  expected 
to  fluctuate  because  of  the  turbulent  eddies  causing 
changes  in  velocity  near  the  wall.  As  a result  of  these 
predicted  effects  at  the  onset  of  turbulence,  an  experi 
ment  was  constructed  to  attempt  to  determine  the  lower 
critical  Reynolds  number  using  streaming  potential  measi 


2.  Laboratory 

The  laboratory  was  especially  designed  and 
constructed  for  the  present  experiments . This  was 
necessary  since  the  work  is  particularly  sensitive  to 
temperature  and  humidity.  The  strong  temperature 
dependence  is  due  to  the  viscosity  and  density  of  the 
fluid  and  even  though  extremely  precise  temperature 


control  was  achieved,  the  changes  In  viscosity  are  still 
the  overriding  errors.  The  humidity  of  the  room  comes 
into  play  in  the  electrical  measurements . The  currents 
involved  are  of  the  order  of  10-8  to  10~9  amps  with  a 
source  Impedance  of  the  order  of  50  meg  ohms . Hence , 
surface  leakage  could  be  a problem. 

The  laboratory  is  a 10  x 16  foot  room  Inside 
an  air-conditioned  building.  The  entire  room  is  vapor 
sealed  with  10  mil  polyethylene.  The  walls  and  floor 
are  insulated  with  3 inches  and  the  ceiling  with  4 inches 
of  Styrofoam.  The  air-conditioning  controls  are  of  the 
potentiometer  type  which  control  a servo-driven  back 
pressure  regulator  installed  in  the  compressor  suction 
line . This  provides  control  of  the  room  temperature  to 
± 0 . 7 5°F  over  the  range  40-70°F.  The  ± 0.75°F  is  a 
continuous  variation  with  approximately  a 4 minute  period. 
Since  the  quantity  of  fluid  in  the  experiment  is  5 liters 
of  water,  the  thermal  inertia  is  large  resulting  in  only 
small  drifts  in  temperature — 0(0. 1°C) — over  a period  of 

The  air  drying  is  accomplished  by  use  of  a 
compression-cooling-expansion  cycle.  This  enables  the 
room  to  be  held  below  20%  R.H.  during  operation. 
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The  apparatus  used  is  schematically  shown  in 
Fig.  2.  The  reservoirs  R^  and  R2  are  5 liter  boiling 
flasks  fitted  with  ground  glass  taper  joints,  and  Jg. 
Tubes  L1  and  are  6 mm  Pyrex  with  taper  Joints  at 
the  top  for  inserting  electrodes.  A ^ and  are  8 mm 
tubes  terminating  with  the  female  portion  of  ground 
glass  ball  joints.  The  tube  T is  a precision  bore 
capillary  tube  0.0242  i .0003  Inches  in  diameter,  4.8 
Inches  long.  The  ends  are  flared  into  8 mm  tubes  which 
end  in  the  male  portion  of  the  ball  joints.  In  operation, 
the  liquid  is  forced  up  L^,  over  a glass  trip  at  B 
(see  Fig.  3)  and  through  T by  applying  gas  pressure 
through  Vj.  The  levels  in  the  flasks  are  determined  by 
connecting  either  or  S2  to  a height  gage  (see  Fig.  4). 
The  level  can  be  determined  within  ± 0.002  inch. 

The  electrodes,  and  Eg,  were  Ag-AgCl 
constructed  on  a spiral  of  24  AffG  platinum.  The  Ag-AgCl 
coating  was  made  with  a mixture  of  90%  silver  oxide  and 
10%  silver  chlorate  formed  into  a paste  and  painted  on  the 
platinum  wire.  The  electrode  was  then  baked  at  500°C  for 


fifteen  minutes. 


Fig.  2 Schematic  sketch  of  experimental 


apparatus. 


Photograph  of  the  glass  trip 
of  the  precision 


Fig.  4 Height  gage  nanometer  used 
determining  the  water  level 


The  gas  pressure  was  supplied  by  commercially 
bottled  nitrogen.  A sketch  of  the  pressure  system  is 
shown  in  Fig.  5.  The  gas  was  reduced  in  pressure  by  a 
standard  pressure  regulator  obtained  from  the  local  gas 
company.  It  was  found  in  the  beginning  that  such  a 
regulator  could  not  control  the  pressure  to  the  desired 
accuracy  so  the  needle  valve  and  the  bleed  valve  were 
installed.  The  bleed  valve  was  designed  for  this 
application  and  is  shown  in  Fig.  6.  This  valve  consists 
simply  of  a 1 mm  capillary  tube  with  provisions  for  a 
26  gage  wire  to  slide  in  and  out  of  the  bore.  The 
pressure  drop  through  the  tube  depends  on  the  position  of 
the  wire.  This  setup  worked  very  well,  providing  a means 
of  controlling  the  pressure  on  the  flasks  very  precisely. 

The  manometer  was  a simple  mercury  manometer  used  in 
combination  with  a microscope  and  a height  gage  (see  Fig.  7) 
With. this  optical  setup,  a change  in  pressure  of  0.001 
inch  of  mercury  could  be  detected.  This  combined  with 
the  precision  bleed  valve  allowed  the  pressure  to  be 
controlled  easily  to  within  ± 0.005  inch  Hg.  The  two 
three-way  stopcocks  in  the  supply  line  provided  a 
convenient  means  of  either  applying  pressure  to  the  flasks 


venting 


Fig.  5 A sketch  of  the  pressure  system. 


Manometer  and  optical  system  used 
for  pressure  measurement. 


difficult  problems 


experiment  was  the  cleaning  of  the  experimental 
components . This  included  cleaning  the  glassware  and 
electrodes  and  maintaining  rather  rigid  purity  controls. 
All  this  was  desirable  since  contaminants  could  cause 
inconsistent  results. 

The  following  procedure  was  used  for  the 
glassware.  The  glassware  was  soaked  in  chromic  acid 
for  two  hours  and  then  rinsed  in  190°P  conductivity 
water.  While  the  glass  was  still  hot,  chromic  acid  was 
added  again  and  allowed  to  stand  for  another  hour.  It 
was  then  rinsed  four  times  in  hot  conductivity  water  and 
steamed  in  the  same  quality  water  for  four  hours.  If 
the  surface  of  the  glass  did  not  wet  uniformly  when 
steaming  began,  it  was  treated  again  with  alternate  doses 
of  chromic  acid  and  hot  water  until  the  wetting  was 
uniform.  Steaming  was  then  continued  for  four  hours. 

The  electrodes  were  cleaned  in  a 10%  solution 
of  nitric  acid  for  two  hours  and  leached  in  conductivity 

The  conductivity  water  was  obtained  from  a 
commercially  built  still.  The  still  was  adjusted  so  the 


condensate  came  out  at  approximately  190°F.  This  was 
collected  In  a 9 liter  Pyrex  serum  bottle  with  a ground 
glass  stopper.  Before  being  used,  the  bottle  was 
cleaned  as  described  above.  During  collection,  nitrogen 
was  bubbled  through  the  condensate.  This  helped  to 
remove  any  remaining  carbon  dioxide  and  maintained  a 
nitrogen  atmosphere  in  the  bottle^  preventing  reabsorption 
of  carbon  dioxide  and  ammonia  from  the  air.  Any  water 
with  a specific  resistance  lower  than  4.0  x 10®  ohm-cm 
was  rejected. 


The  solution  used  in  the  experiments  w 
0.001  M KC1,  which  resulted  in  a double  layer  t 
of  approximately  10"®  cm.  This  was  made  using  conduc- 
tivity water  and  reagent  grade  KC1  which  had  been 
recrystallized  in  the  same  quality  water.  The  mixing  of 
this  solution  was  timed  so  that  it  could  be  placed  in  the 
flasks  as  soon  as  the  steaming  process  had  ended.  The 
entire  apparatus  was  then  sealed  immediately.  The 
sealing  was  accomplished  using  red  sealing  wax;  the  wax  w 
never  allowed  to  come  into  contact  with  the  solution. 

The  flasks  were  never  reopened  until  the  completion  of 
the  experiments . The  procedure  seemed  to  be  successful 
in  that  the  streaming  potentials  were  quite  consistent 


Electrical  Measurements 


The  electrical  measurements  were  afflicted  with 
the  usual  problems  found  in  electrochemistry,  i.e.,  low 
currents,  0(10”  amps),  and  high  source  impedance, 

0(50  meg  ohms) . Thus  the  measuring  device  selected  had 
very  low  current  drain  and  the  entire  experiment  was 
shielded  from  stray  electrical  noise.  A schematic 
diagram  for  the  measurements  is  shown  in  Fig.  8.  The 
experiment  was  enclosed  in  a screen  room  which  was 
grounded  at  G by  a 20  foot  ground  rod.  A General  Radio 
Model  number  1230- A Electrometer,  having  an  input 
impedance  of  10*1  ohms,  was  used  as  a preamplifier.  The 
electrodes  were  connected  to  the  electrometer  by  a coaxial 
conductor  passing  through  the  wall  of  the  screen  room. 

The  output  of  the  electrometer  was  read  in  two  ways 
depending  on  the  experiment . If  the  average  value  of  the 
streaming  potential  was  desired,  the  output  was  connected 
directly  to  a Dymec  2401A  integrating  digital  voltmeter. 

If  the  fluctuations  of  the  streaming  potential  were  to  be 
studied,  the  output  voltage  of  the  electrometer  was  biased 
by  a d.c.  voltage  so  the  resulting  signal  could  be 
recorded  by  a 10  inch  Varian  strip  chart  recorder  operating 
on  a full  scale  range  of  5 rav.  This  arrangement  proved  to 
be  quite  satisfactory.  The  results  are  shown  in  Chapter  III. 
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Electrical  diagram  for 
and  shielding. 
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6.  Experimental  Procedures 

The  final  procedure  used  to  determine  the 
critical  Reynolds  number  was  to  measure  the  pressure 
where  a particular  type  of  fluctuation  (to  be  discussed 
in  Chapter  III)  in  the  streaming  potential  was  found. 
This  pressure  was  then  used  to  find  the  corresponding 
Reynolds  number  from  calibration  curves  obtained  from 
mass  flow  measurements.  The  calibration  curve  is 
shown  in  Fig.  9.  This  curve  is  a linear  least  squares 


the  indicated 


points . 


III.  RESULTS 


The  first  measurements  taken  were  the  average 
streaming  potential  as  a function  of  the  pressure 

Pd ’ 

differential  causing  the  flow.  A typical  plot  of  ja-jX 
versus  E is  shown  in  Pig.  10  where 
P is  the  pressure, 
d is  the  diameter  of  the  tube, 
ja.  is  the  absolute  viscosity, 

1)  is  the  kinematic  viscosity, 

J.  is  the  length  of  the  tube, 
and  E is  the  streaming  potential. 

It  is  seen  that  there  seems  to  be  a slight  break  in  the 
curve  in  the  region  = 10.30  This  is  emphasized 

by  the  two  straight  lines  drawn  through  the  data  points. 
However,  the  break  is  so  slight,  it  is  difficult  to 
select  accurately  the  pressure  where  this  occurs.  To 
avoid  this  problem,  the  data  were  taken  at  pressure 
Intervals  of  0.1  inches  of  Hg  (Instead  of  0.5  inches  as 
used  in  Pig.  10)  and  then  differentiated  one  time.  In 


■ 


this  way,  it  was  thought  the  break  could  be  detected  with 
a great  deal  more  precision.  The  method  was  unsuccessful. 

It  turned  out  that  in  the  region  of  interest  the  sample 
period  of  the  integrating  voltmeter  was  too  short  (1  sec.) 
to  give  a steady  reading  to  the  desired  accuracy.  The 
readings  would  oscillate  roughly  - 0.1  mv.  Since  the 
change  in  streaming  potential  was  of  the  order  of  10  mv/ 
inch  Hg,  this  represented  an  error  of  ± 10%  in  the 
measured  slope  of  the  curve  for  pressure  increments  of 
0.1  inch.  The  error  completely  obliterated  the  change 
in  slope  at  the  critical  Reynolds  number. 

A decision  was  made  to  use  this  troublesome 
noise  rather  than  try  to  eliminate  its  effects.  To  do 
this,  it  was  necessary  to  look  at  a small  noise  component 
which  had  a characteristic  peak-to-peak  amplitude  of  0.25  mv 
and  which  was  additionally  superimposed  on  a d.c.  voltage 
of  the  order  of  125  mv.  This  was  accomplished,  as  indicated 
in  Chapter  II,  by  biasing  the  signal  with  an  opposing  d.c. 
voltage  so  that  the  noise  could  be  recorded  on  a 10  inch 
strip  chart  recorder  with  a 5 mv  full  scale  deflection. 
Typical  records  are  shown  in  Figs.  11,  12,  13  and  14  for  a 
temperature  of  19.76°C.  Traces  1 and  2 of  Fig.  11  show 


Fig.  11  Typical  traces  for  low  flow  rates. 
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the  noise  level  for  pressures  of  12.83  and  14.33  inches 
of  Hg.  No  change  in  the  traces  can  be  seen  from  zero 
flow  up  to  a region  just  below  transition.  Traces  3 and 
4 of  Fig.  12  show  a gradual  increase  in  noise  up  to 
transition. 

First  thoughts  were  that  the  beginning  of  this 
buildup  represented  instability  in  the  flow.  However, 
the  fact  that  it  increased  gradually  could  not  be  explained 
in  the  light  of  the  results  of  earlier  experimenters  who 
observed  the  onset  of  turbulence  to  be  quite  rapid  and 
violent.  Thus,  the  search  was  continued  at  higher  and 
higher  pressures  to  see  if  there  was  anything  to  Indicate 
a sudden  change  in  affairs  within  the  flow.  Indeed,  such 
an  indication  was  found.  This  is  shown  by  traces  5 and  6 
in  Fig'.  13.  The  sharp  spikes  rising  well  above  the  back- 
ground appeared  very  suddenly  and  became  more  numerous  as 
the  pressure  was  Increased.  A very  curious  thing  to  be 
noted  is  that  spikes  always  represented  an  increase  in 
voltage  when  they  first  appeared.  As  the  pressure  is 
Increased  to  a sufficiently  high  value,  the  spikes  become 
positive  and  negative  and  completely  dominate  the  trace. 


Typical  traces  showing  the  buildup 
of  noise  as  transition  is  approached. 


Fig.  13 


noise  and  the  spikes  caused  by  instability 
in  the  flow. 


Fig.  14  Typical  trace  for  flow 
critical  state. 


Having  completed  this  study,  it  was  reasoned 
that  the  initial  noise,  shown  in  traces  3 and  4,  was 
caused  by  disturbances  from  the  trip  at  the  mouth  of  the 
tube  which  had  not  been  completely  damped  before  reaching 
the  wall.  The  velocity  at  which  the  high  peaks  first 
began  was  taken  as  the  critical  velocity.  Using  this 
criterion,  32  measurements  were  made  at  various  tempera- 
tures to  determine  the  critical  Reynolds  number . The 
results  are  tabulated  in  Table  1.  The  pressure  is 
included  in  the  table  in  its  dimensional  form  to  indicate 
the  sensitivity  of  the  results  to  temperature  variations. 
Notice  that  a temperature  change  from  20.54°C  to  19.50°C 
required  a change  in  pressure  from  13.68  to  16.63  Inches 
of  Hg  in  order  to  reach  transition.  The  mean  value  of 
the  Reynolds  number  was  found  to  be  1907  with  a mean 
error  of  - 3.27.  The  standard  deviation  of  the  measure- 


1907  ± 3 is  in  excellent  agreement  with  1900  found  by 
Barnes  and  Coker  (loc.  cit.)  and  with  the  recent  work  of 
Lessen,  Fox,  Bhat  and  Lin^^  where  a range  from  1400  to 


reported. 


1949 


15.50 

16.05 


1915 


1890 

1894 
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